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Abstract 

Considering charged fluid spheres as anisotropic sources and the dif- 
fusion limit as the transport mechanism, we suppose that the inner space- 
time admits self-similarity. Matching the interior solution with the Reissner- 
Nordstrom-Vaidya exterior one, we find an extremely compact and oscil- 
latory final state with a redistribution of the electric charge function and 
non zero pressure profiles. 
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1 Introduction 

It is well known that astrophysical objects are not significantly charged, For 
this reason charged bodies have limited interest (for an historical overview see 
Section I in Ref. pQ). Observed stars, like the Sun or a neutron star, cannot 
support a great amount of charge. This is so basically because the particles that 
compose a star have a huge charge to mass ratio, as is the case for a proton or 
an electron. For cold stars the electric charge is about 100 Coulombs per solar 
mass [2]. Nevertheless, in compact stars the equilibrium conditions are different, 
allowing some more net charge [3]. For highly relativistic stars near to form a 
black hole, the huge gravitational pull can be balanced by huge amounts of net 
charge. On this kind of configurations we are concerned in this paper and have 
been considered by other authors PQ,[3]-[7]. 

To include the electric charge a number of authors make additional assump- 
tions such as an equation of state or a relationship between metric variables [8]- 
PTj . Bonnor and Wickramasuriya [T^] have studied electrically charged matter, 
with electrostatic repulsion balancing the gravitational attraction. Most works 
have been done under static conditions, including the ones by Ivanov |13j . who 
exhaustively surveyed static charged perfect fluid spheres in general relativity, 
and Ray et al. p], who studied the effect of electric charge in compact stars 
and its consequences on the gravitational collapse. Bekenstein [3] found that 
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for highly relativistic stars, whose radius is on the verge of forming an event 
horizon, the large gravitational pull can be balanced by large amounts of net 
charge. 

In this paper we study an electrically charged matter distribution as an 
anisotropic fluid. It is well known that different energy-momentum tensors 
can lead to the same space-time |14j . For instance, under spherical symmetry 
viscosity can be considered as a special case of anisotropy [15] . To illustrate our 
approach to obtain dynamical models we consider the diffusion approximation 
as the transport mechanism, and a self-similar space-time for the inner region. 
Our treatment allows different scenarios for the gravitational collapse, including 
the reported some years ago 16 . In particular we are now interested in non 
stationary initial conditions to explore the fate of the gravitational collapse. The 
set problem has some physical and geometrical features that we review briefly, 
these are: anisotropy, dissipation and self-similarity. 

For certain density ranges, local anisotropy in pressure can be physically 
justified in self-gravitating systems, since different kinds of physical phenomena 
may take place, giving rise to local anisotropy and in turn relaxing the upper 
limits imposed on the maximum value of the surface gravitational potential [17] . 
The influence of local anisotropy in general relativity has been studied mostly 
under static conditions (see [8] and references therein). Herrera et al. [TB] have 
reported a general study for spherically symmetric dissipative anisotropic fluids 
with emphasis on the relationship among the Weyl tensor, the shear tensor, 
the anisotropy of the pressure and the density inhomogeneity. In another con- 
text, a scalar field with non-zero spatial gradient is an example of a physical 
system where the pressure is clearly anisotropic. Boson stars, hypothetical self- 
gravitating compact objects resulting from the coupling of a complex scalar field 
to gravity, are systems where anisotropic pressure occurs naturally. Similarly, 
the energy-momentum tensor of both electromagnetic and fermionic fields are 
anisotropic 19J. We want to stress that in order to have isotropy we need an 
extra assumption on the behavior of the fields or on the fluid describing interiors. 

The proposed dissipative process is suggested by the fact that they provide 
the only plausible mechanism to carry away the bulk of binding energy, leading 
to a neutron star or black hole [5D]. Furthermore, in cores of densities close to 
10 12 g / cm? the mean free path for neutrinos becomes small enough to justify 
the use of the diffusion approximation [211 E2]- Here we do not discuss the 
temperature distribution during diffusion; we consider this transport mechanism 
because it provides an efficient way to radiate energy. Additional studies are 
necessary to explore the consequences of incorporating a hyperbolic theory of 
dissipation to overcome the difficulties inherent to parabolic theories. 

Few exact solutions to the Einstein-Maxwell equations are relevant to grav- 
itational collapse. For this reason, new collapse solutions are very useful, even if 
they are simplified ones [33] . It is well known that the field equations admit ho- 
mothetic motion [25j-|27j. Applications of self-similarity range from modeling 
black holes to producing counterexamples to the cosmic censorship conjecture 
[28]-[35]. It is well established that in the critical gravitational collapse of an 
scalar field the space-time can be self-similar [36]-[38]. We can expect on phys- 
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ical grounds that this scenario remains similar for charged matter [3T)] . 

In this work the Darmois-Lichnerowicz junction conditions at the surface of 
the distribution are satisfied. We match the self-similar interior solution with 
the Reissner-Nordstrom-Vaidya exterior one. We obtain a model that resembles 
the Seidel-Suen solitonic star formed by a real and massive scalar field [40 . 

Section [2] contains the field equations, written in such form that the electri- 
cally charged fluid can be considered anisotropic. In this section we also discuss 
the junction conditions to match the dissipative interior space-time with the 
exterior space-time, a Reissner-Nordstrom-Vaidya one, and write the gener- 
alized Tolman-Oppcnhcimcr-Volkoff equation. Section [3] contains a review of 
self-similarity for the sake of completeness. Then we propose a simple solution 
to include electric charge by means of the active mass. We show two example 
solutions in section [4] and finally discuss our results in section [5] 



2 The field equations and Junction conditions 

Let us consider a non static distribution of matter which is spherically symmet- 
ric and consists of charged fluid of energy density p, pressure p, electric energy 
density e and radiation energy flux q diffusing in the radial direction, as mea- 
sured by a local Minkowskian observer comoving with the fluid (with velocity 
— u>). In radiation coordinates |H] the metric takes the form 

ds 2 = e 213 ^— du 2 + 2du dr^j - r 2 (d6 2 + sin 2 6 d(j) 2 ^j , (1) 

where j3 and V are functions of u and r. As it is well known for the spherical 
symmetry F ur = —F ru are the only non vanishing electromagnetic field tensor 
components. Now, defining the function C(u,r) by the relation 

F ur = £ e -2/3 ( 2 ) 

the inhomogeneous Maxwell equations become 

C r = 4irr 2 J u e 213 , (3) 

and 

C tU = -47rr 2 J r e 2 ^ (4) 

where the comma subscript represents partial derivative with respect to the in- 
dicated coordinate; J u and J r are the temporal and radial components of the 
electric current four vector, respectively. Thus, the function C(u,r) is natu- 
rally interpreted as the charge within the radius r at the time u. Thus, the 
inhomogeneous Maxwell field equations entail the conservation of charge 

tf M C,M = 0, (5) 
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where the four-velocity is given by 
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We can write the conservation equation in a more suitable form, that is 

C,u + %-C, r = 0, (7) 
du 

where 



dr 
du 



V OJ 



is the matter velocity. It is clear that inside a sphere comoving with the fluid 
the charge is conserved. 

Introducing the mass function by means of 



M = \ (r - Ve^) 



we can write the Einstein field equations as follows 
p + p r u> 2 2ujq 
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with p = p + e, p r =p — e, p t =p + e and the electric energy density e = E 2 /8tt, 
where E = C/r 2 is the local electric field intensity. Observe that we have 
introduced the mass function p instead the usual total mass m = (r— Ve~ 2,3 )/2+ 
C 2 /2r. 

If we define the degree of local anisotropy induced by charge as A = pt—p r = 
2e, the electric charge determines such a degree at any point. 



4 



It can be shown that the Tolman-Oppenheimer-Volkoff equation reduces to 



dp , P + P 
dr 1 — 2/i/r 
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(14) 



This generalized equation is the same as for an anisotropic fluid [43] . 

The most general junction conditions for our system can be written as [TC] 



Ma 



Crp 

2a ' 



/3a = 0, 



and 
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(15) 
(16) 
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2a 4a 3 

where a subscript a indicates that the quantity is evaluated at the surface, rh a 
and Ct are the total mass and the total charge, respectively. Remarkably, the 
last equation is equivalent to 

Pa = q a - (18) 

We have obtained a system of equations which are equivalent to the anisotropic 
matter case (see [5] and references therein). In the next section we will see how 
this approach lead us to an extra simplification. 

Up to this point all the written equations are general. We can try to solve the 
system using numerical techniques or a seminumerical approach to gain some 
insight. The last method is worked out in the next section. 



3 The self— similar space— time: A simple solu- 
tion 

We shall assume that the spherical distribution admits a one-parameter group 
of homothetic motion. A homothetic vector field on the manifold is one that 
satisfies £^g =2ng on a local chart, where n is a constant on the manifold and 
£ denotes the Lie derivative operator. If n ^ we have a proper homothetic 
vector field and it can always be scaled so as to have n = 1; if n = then £ is a 
Killing vector on the manifold. So, for a constant rescaling, £ satisfies 

= 2g (19) 

and has the form 

£ = A(u,r)d u + \{u,r)d r . (20) 
The homothetic equations reduce to 

m = o, (2i) 
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£(w) = 0, (22) 

where A = r, A = A(u), x = A*/ 7 " an( i v = Ae 2/3 /r. Therefore, x — x(C) an d 
v = v(Q are solutions if the self-similar variable is defined as 

( = re-f du/A . (23) 

We have reported [HI HI] a simple homothetic solution which can be written 
for the present case as 

li = tia{u){r/a) k+1 (24) 

and 

e 2 ? = (r/a) l+1 (25) 

with k and I constants. It can be shown that for a given k, I must be a root of 
a complicated seventh degree polynomial in order to satisfy the additional sym- 
metry equations. It is obvious that these solutions satisfy the continuity of the 
first fundamental form, and to satisfy the continuity of the second fundamental 
form, the radial velocity at the surface must be 

2(Z + l)(l-2/i a / a ) 
2(fc+l)/i a /a-C 2 /a 2 



3.1 Surface equations 

Now, to completely determine the metric we have two differential equations at 
the surface, the equation ([8]) and the field equation (fT0|) which evaluated at the 
surface become respectively 

d = (l-2 Ma /a)-^- (27) 

J- — LL> a 

and 

C 2 

Aa = -0(1 - 2/i a /o) + O^, (28) 

where a dot over a variable denotes the derivative with respect to time. The 
quantity Q is defined as 

Q=^a\ a (\±^\, (29) 

V 1 - 0J a J 

and can be explicitly determined using equation IT4]) evaluated at the surface, 
resulting in 

Q = {(1 - 2^ Q /a)[2(; + l) 2 + k(k + 1 + 3(1 + 1))] - 2(C T /a) 2 
-k(k + 1 + 3(/ + 1))}{A + 1 - (1 - 2 Ala /a)(fc + I + 2) - (C T /a) 2 }/ 

[4(1 - 2fi a /a)[l + l)a 2 ] . (30) 

The surface equations (|27|) and (|28]) together with equations (|26|) and p0|) 
can be integrated numerically (using Runge-Kutta for instance) for initial con- 
ditions a(0) and (J. a (0) and some k and Ct- The only restrictions on these 
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conditions and parameters come from the need to get a physically acceptable 
model. 

For k — (which implies I = 0) we obtain the same results of reference |16j . 
As we have mention before, we gain some simplification with the anisotropic 
approach of this paper, that is, one of the similarity equations (Eq. (23) in Ref. 
[IB] ) was decoupled, leading us to less restrictive models. 

3.2 Integrating the conservation equation 

Once the surface equations are integrated we have to integrate the conservation 
equation to obtain all the physical variables inside the source. To do that 
we use comoving spatial markers x = r/a. Thus the conservation equation can 
be written as 

C u = j~C ' x (31) 
du 

which is a wave-like equation that can be integrated numerically using the Lax 
method (with the appropriate Courant-Friedrichs-Levy (CFL) condition). The 
evolution of the conservation equation is restricted by the surface evolution and 
was implemented as follows 

cr - 5 + q-o + £ (£)' - q-o ■ m 

The superscript n indicates the hypersurface u = nSu and the subscript j the 
spatial position for a comoving observer at x = jdx. We typically used Su = 
10~ 2 with a CFL condition Su = 25x. In order to integrate the conservation 
equation we need to specify a boundary condition and an initial condition. For 
this work we have used the boundary condition 

C(x = 0,u)=0, (33) 

with an initial condition 

C{x,u = 0) = Ctx v , (34) 

where the power V allows us to test the sensitivity to the initial conditions of 
our results. 



4 Modeling 

4.1 Evolving towards dust: k = 

As a test of our approach to dealing with electric charge as anisotropy is to 
reproduce the results reported in [16] , For this scenario k = I = 0, a(0) = 3.5, 
A*a(0) = 1, V = 1 and Ct — 0.09. Even when the initial mass m a (0) is not 
equal to one, we obtain the same results, that is, the collapse is halted and the 
boundary oscillates when the total electric charge asymptotically approximates 
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the total mass m a and the radius equals asymptotically to twice the total electric 
charge. The ratio P/p as a function of time is the same for any region. Also, as 
the distribution evolves, becomes dust asymptotically with damped oscillations 
and the pressure at the surface coincides with the heat flow. The matter velocity 
profiles at all regions have damped oscillations that evolve to the static regime. 
A striking result not reported before is the ability to obtain any interior profile 
for matter velocity or electric charge function using the simple rule 

!F{u, x) = xJ-(u, 1), (35) 

where T could be either dr/du or C, that is, once the surface profile is de- 
termined we can go without additional work inside the distribution. Another 
interesting result is that we can get the interior profiles for pressure and energy 
density by 

g( u ,x) = g(u,i)/x 2 , (36) 

where Q could be p or p. 

4.2 Redistribution of the electric charge: k = 0.1 

Now that we have successfully tested our new approach, we explore another 
physical scenario, one with redistribution of the electric charge. For this case 
we choose k — 0.1, a(0) = 3.5, p a {0) = 1 and Ct = 0.09 to obtain (among 
other roots) I — 0.741632. For V = 1 (power of the initial condition to integrate 
the conservation equation) the results are shown in figures HHS The system 
collapses and losses mass until it reaches a stable oscillatory state. When the 
distribution is in oscillatory regime the mass behaves in the same way, that is, 
oscillates between two extremal values. This oscillatory behavior extends to all 
physical variables at all pieces of the material. The interior redistribution of 
electric charge responds to the surface dynamics. Electric charge (repulsively) 
let the distribution expand with absorption of energy up to some maximum 
point in which gravitation compensates electric repulsion. Therefore the collapse 
undergoes with emission of energy, up to some minimum point in which electric 
repulsion is again dominant and so on indefinitely. This model seems to be 
independent of the initial conditions, as shown in figure [7J The dependence of 
the radius on the total electric charge is shown in figure [8] Analysis of this 
last result indicates scale invariance with respect to the total electric charge. In 
fact, once the stationary state (oscillatory) is reached, the amplitude and period 
depend on the total electric charge. Thus, if we know the oscillatory regime for 
some specific value of electric charge, we can reproduce the oscillatory regime 
for any other, non zero, value of electric charge. In the limit of zero charge 
the distribution shrink completely, forming a very small region with very high 
energy density and pressure. This last result was found in another context (see 
section IV. B in Ref. [33]). 

The energetic at the boundary surface in conjunction with the pull of gravity 
and push of electric charge, explain the dependence between the amplitude 
and period with the total electric charge. These features are analogous to the 
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reported by Brady et al. [12] for a real and massive scalar field. In this context 
the mass destroys the scale invariance of the Einstein-scalar field equations, but 
after the nearly critical phase the oscillation period depends on the inverse of 
scalar field mass, i.e. the Compton wavelength. In our case the situation is 
analogous respect to the electric total charge. Parenthetically, the Seidel-Suen 
oscillaton [10] is the stable solution far from the critical behavior of a massive 
and real scalar field. Our solution behaves like that stable configuration from 
quite simple considerations. 

5 Discussion 

In this paper we have shown that electrically charged matter can be considered 
as anisotropic matter. We exemplify the approach obtaining a dynamical model 
under the diffusion approximation and self-similarity within the source. The 
showed example is representative of many others varying only the initial condi- 
tions, k and Ct, and imposing the following physical conditions: — 1 < lu < 1, 
p < p and p > 0. 

Some general considerations concerning our results are listed next, 

• Extremely compact distribution: The collapse is halted when the 
total electric charge is about the total mass (in geometrical units). If the 
initial radius of the electrically charged distribution is 5.2 km, the radius 
of the final stable source is » 266 m. Also the final pressure (if not zero) 
is about 1.5 x 10 39 Pascal and the energy density 1.8 x 10 19 g/cm 3 . 

• High electric field: The electric charge used in our models can be 
as high as 1.5 x 10 19 Coulomb. Thus, the electric field for the stable 
configuration is about 2.2 x 10 19 V/m P2Q]. The sign of net charge does 
not affect our results. 

• High powered: If the initial total mass is 1M the final is w 0.1M Q , 
that is, almost 90% of the the mass is radiated in 0.15 ms. 

A more general solutions to the homothctic equations ([2"T]) and (j2"2"|) could 
be necessary to explore how general our results are. However, the results found 
using the anisotropic approach are consistent with previously reported results 
when the static regime is the final state [45]-[49]. If the solution reaches sta- 
tionarity it oscillates. For this case we did not investigate the stability of the 
system, but independence of initial conditions is apparent. It is clear that the 
presence of some finite electric charge (or anisotropy) avoids singularity forma- 
tion. Finally, we would like to conclude by posing a question. If the spherically 
symmetric source has an upper limit for the total net electric charge that carries 
[TJ [SU] , is it relaxed by the transport mechanism or by the additional homothctic 
symmetry? Work in this direction is in progress. 
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Figure 1: Total electric charge (short dashed line) and the evolution of radius 
a (solid line) and the total mass m a (dashed line). 
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Figure 2: Energy density as a function of time for different pieces of the dis- 
tribution: x — 0.25 (short dashed line); x — 0.5 (large dashed line); x = 0.75 
(solid line) and x — 1 (dotted line). 
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Figure 3: Pressure as a function of time for different pieces of the distribution: 
x = 0.25 (short dashed line); x — 0.5 (large dashed line); x — 0.75 (solid line) 
and x = 1 (dotted line). 
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Figure 4: Heat flux as a function of time for different pieces of the distribution: 
x — 0.25 (short dashed line); x — 0.5 (large dashed line); x — 0.75 (solid line) 
and x = 1 (dotted line). 
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Figure 5: Matter Velocity as a function of time for different pieces of the dis- 
tribution: x — 0.25 (short dashed line); x — 0.5 (large dashed line); x = 0.75 
(solid line) and x = 1 (dotted line). 
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Figure 6: Electric charge function versus time for different pieces of the material: 
x = 0.25 (solid line); x = 0.5 (large dashed line); x = 0.75 (short dashed line). 
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Figure 7: Electric charge function versus time for x = 0.75 and different values 
of power V in the initial value of C to integrate the conservation equation. 
Curves correspond to: V = 1 (continuous line); V = 2 (large dashed line); 

= 3 (short dashed line) and V = 4 (dotted line). Many values of a(0) were 
used with essentially the same results. 
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Figure 8: Radius as a function time for CV = 0.00 (solid line); Ct — 0.09 (large 
dashed line); Ct — 0.18 (short dashed line); Ct — 0.36 (dotted line). 
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